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Abstract 

In continuous-time wavelet analysis, most wavelet present some kind of symmetry. Based on the Fourier and Hartley transform 
kernels, a new wavelet multiresolution analysis is proposed. This approach is based on a pair of orthogonal wavelet functions and is 
named as the Fourier-Like and Hartley-Like wavelet analysis. A Hilbert transform analysis on the wavelet theory is also included. 
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1 Introduction 

The idea of comparing Fourier analysis with wavelet decompo¬ 
sitions is the starting point for introducing an analysis based on a 
couple of orthogonal wavelet functions, one with even symmetry 
and other with odd symmetry. This approach is here presented 
and named as the Fourier-Like and Hartley-Like wavelet analy¬ 
sis. In the standard Fourier analysis, a signal x(t) is simultane¬ 
ously analyzed by even and odd (in quadrature) functions, being 
represented by: 

x(t) = d.c. term + cosine terms + sine terms. 

In the standard wavelet multiresolution analysis (WMRA) (TJ, 
x(t) may be represented by: 

x(t) = (j) term+ i jf terms, 

where “0 term” account for the analysis of x(t) with a scaling 
function (j>{t) and the jf terms” represent those ones derived 
from scaled versions of a mother wavelet function i f{t). 

Comparing the WMRA and the Fourier analysis equations, 
the scaling coefficients, “0 term, play a role that corresponds 
to the d.c. term of the Fourier series; the wavelet coefficients, 
“i f/ terms”, can be viewed as the harmonic components of the 
Fourier series, since the harmonics are scaled versions of the in¬ 
finite Fourier kernel. 

Hence, it seems natural to associate to each wavelet with 
even symmetry another one of odd symmetry, and vice-versa. 
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The Hilbert transform can naturally be invoked to derive the 
in quadrature version of a symmetric (or an anti-symmetric) 
wavelet. Hence, x(t) may be represented on a new continuous¬ 
time WMRA by: 


x(t) = 0 term + i jf terms + orthogonal of y/ terms. 

Accordingly, new wavelet functions that look like the Fourier 
and Hartley transform kernels are invoked. So, by analogy to 
Fourier and Hartley transforms, the cosine and sine kernel is re¬ 
placed by yr and the Hilbert transform of i jf in this new concept 
of wavelet analysis. In order to allow further investigation on the 
Fourier-Like and Hartley-Like wavelet analysis, a brief review of 
the Hilbert transform is presented, as well the results of applying 
some wavelet properties to the Hilbert transform of wavelets. 


2 The Hilbert Transform 


The Hilbert transform of a function g(t) is defined by j2j: 

^{s}(0=P-'A/ r~j: dT ’ 

71 . 7-00 t — T 

where p.v. is the Cauchy principal value of the integral. After 
a change of variable, the Hilbert transform can be written as a 
convolution: 

^{8}(t) = l_*g(t). ( 1 ) 

7T 

By taking the Fourier transform of ([!]), we have: 


where & is the Fourier transform operator and sgn(-) is the 
signum function defined by: 


sgn(ffl) 


1 , ©> 0 , 

< 0 , co = 0 , 

— 1, (0<O. 


( 3 ) 
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It follows from ([2])-([3]) that the Hilbert transform of a func¬ 
tion imposes a null at co = 0 and a — Ti/l phase delay on the 
frequency response of that function. Other interesting properties 
of the Hilbert transform are ||2): 

• A function and its Hilbert transform are orthogonal over the 
infinite interval; 


Proof: The nth moment of \j/(t) is defined by 0 MnM)\ = 
S-oo tn • As i j/(t) has N vanishing moments, then 

Mb[w(t)\ = 0, for n = 0,1 ,...,N — 1. In the frequency domain, 

the moments of \ff(t) are expressed by ij^j: M n [\ff(t)\ = ’ 

where the superscript (n) denotes the nth derivative of T'(cd). 
Hence, the nth moment of f}(t) is given by 


• The Hilbert transform of a real function is a real function; 

• The Hilbert transform of an even function is an odd func¬ 
tion, and vice-versa. 

In the framework of wavelets, the Hilbert transform of a symmet¬ 
ric (or anti-symmetric) real wavelet is a real anti-symmetric (or 
symmetric) function. However, it is necessary to verify whether 
the resulting function is also a wavelet function or not. 

3 The Hilbert Transform on The Wavelet Analysis 


M n [Jf{ W }(t)} = 


[—j •sgn(co) • < P(q))]M 

(2 njY 


G)=0 


Consequently, it follows that J4?{\j/}(t) has at least N vanishing 
moments. □ 

In view of Propositions [I] and [2j it follows that J4?{\f}(t) is a 
wavelet with same energy, admissibility coefficient and at least 
same number of null moments than its generating wavelet y(t). 
Figure [T] shows a few continuous-time real wavelets and their 
corresponding Hilbert transforms. 


A function i jf(t) is a mother-wavelet, if and only if, (i) y(t) is 
in the space of finite energy functions L 2 (M), and (ii) 1 j/(t) satis¬ 
fies the admissibility condition [ 11. Some properties are explored 
in the following propositions in view of applying Hilbert trans¬ 
form to wavelets. Let T^co) be the Fourier transform of \//(t), 
be the Hilbert transform of y r(t), (o[\\f(t)] as the en¬ 
ergy of i j/(t), and C[\ff(t)] as the admissibility coefficient of 1 //(t). 

Proposition 1 If \p{t) is a real wavelet, then J4?{\i/}(t) is also 
a real wavelet with same energy and admissibility coefficient of 
the generating wavelet 1 f/(t). 

Proof: If \jf(t) is a real wavelet, then \\f(t) belongs to L 2 (M) and 
satisfies the admissibility condition. Invoking Parseval’s theo¬ 
rem, the energy and admissibility coefficient of J(?{\f}(t) are 
given by: 

= 3-^|-y. S gn(ffl)-»P(ft))| 2 dw, 


4 The Fourier Kernel on The Wavelet Analysis 

It is time to define a new wavelet function that looks like the 
Fourier transform kernel, which can analyze both symmetries of 
an asymmetric signal. The Fourier transform kernel, or Fourier 
kernel, is defined by e^ = cos(t) + j • sin(t). Considering that 
J$?{cos}(t) = — sin(f), then the Fourier kernel can also be writ¬ 
ten as e jt = cos (t) — j J4?{cos}(t). 

This naive observation motivates the definition of Fourier-Like 
wavelets based on a real wavelet and its Hilbert transform. Let 
us define the Fourier-Like wavelet, T{i J/(t)}, by: 

¥(t)} = -j= ■ [V(t) -j-M’iV}(t )}. 

The following proposition demonstrates that T{is also 
a wavelet and that the factor 1 /is imposed so as to guarantee 
that the Fourier kernel holds the same energy and admissibility 
coefficient of its generating wavelet. 

Proposition 3 If\j/(t)isa real wavelet and \f/} (t) its Hilbert 
transform, then Fty/(t ) is also a real wavelet with same energy 
and admissibility coefficient of its generating wavelet\\f{t). 


A simple manipulation yields to 

1 r°° 

c[jrwm = jTJ^d«. 

Applying Parseval’s theorem on the right-side of the energy 
equation, it is straightforward to conclude that d^f{\\f}{t) is 
also in L 2 (M). Moreover, \jf(t) and J^{\j/}{t) have the same 
energy. As i j/(t) G L 2 (M) and 1^(0) | = 0, it promptly fol¬ 
lows that J(?{\j/}(t) also satisfies the admissibility condition: 
C[J^{\\f}(t)] < +oo. Furthermore, \\f(t) and J4?{\\f}(t) have the 
same admissibility coefficient. □ 


Proof: If 1 jf(t) is a real wavelet and J^{\j/}(t) its Hilbert trans¬ 
form, then i j/(t) and Jtf{\\f}(t) belong to L 2 (M) and satisfy the 
admissibility condition. The energy and admissibility coefficient 
of T{y/f7)} are given by: 

1 r°° 

Simple manipulations furnish: 


1 r°° 

#mvm = 2jjv(t)\ 2 +mvm 2 dt, 


cpivK 0 ] 


2 -|'P(<a) 

1 M 


dft), 


co > 0 
co < 0. 


Proposition 2 Let y/(t) be a wavelet with N vanishing moments, 
then y/}(t) has at least N vanishing moments. 
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Figure 1: Continuous-time wavelets and their Hilbert transforms: (a) Morlet; (b) Meyer; (c) Mexican Hat; (d) Gaussian-1; 
(e) Gaussian-2; and (f) Gaussian-3. Scale of the horizontal axis is one unit (time) per division. Scale of the vertical axis is 0.2 
units (amplitude) per division, ranging from — 1 to 1. 


From Proposition [I] it follows that also belongs 

to L 2 (R) and it has the same energy as 1 j/(t). Since \\f(t) G L 2 (R) 
and |*F(0) | = 0, it follows that J{\\f(t)} also has the same admis¬ 
sibility coefficient of i f{t). □ 

The next proposition shows that Fti \f/(t) has same vanishing 
moments than its generating wavelet, y(t). 

Proposition 4 Let \\f(t) be a wavelet with N vanishing moments, 
then 3{y}(t) has [[-Oalso N vanishing moments. 


Proof: From Proposition 2, it follows that the nth moment of 
J{\l/}(t) is given by: 


MninVKt)] 


['P(ffl) — sgn(ffl) •'P(fi))] (n) 

( V 2)»-(-2 -n-j)* 


which can also be written as 


si {g}(0 is a complex signal designed by a real function g(t) and 
its Hilbert transform Jff{g}(t) |3|. In the framework of wavelets, 
an analytic wavelet, s/{\ff}(t), can be defined by: 

**W(0 = {¥}(())■ 

Analytic wavelets have also same energy, admissibility coeffi¬ 
cient and null moments than their generating wavelet, \\f(t). The 
proofs are similar to that of Propositions [3] and [4] In the fre¬ 
quency domain, analytic wavelets are null for co < 0. For co > 0, 
they have the magnitude response of the generating wavelet mul¬ 
tiplied by a scalar factor. 

4.2 Computing the Wavelet Analysis of Asymmet¬ 
rical Real Signals 


M„[TM(0] = M n [\if(t)] + yM n [<^M(0] • 

Then, 3{\\f}(t) has also N null moments. □ 

In the frequency domain, Fourier-Like wavelets are null for 
co > 0. For co < 0, they have the magnitude response of the 
generating wavelet multiplied by a scalar factor. This is a typical 
behavior of analytical signals. 


It may be expected that using Fourier-Like or analytic wavelets, 
even and odd parts of an asymmetric real signal can be better 
analyzed, respectively, by an even wavelet and its Hilbert trans¬ 
form, i.e. an odd wavelet. In both cases, it will be necessary to 
perform a complex wavelet analysis. It is also possible to ana¬ 
lyze both symmetries of a real signal using a real wavelet. In that 
case, the Hartley kernel should be invoked. 


4 1 Analytic Wavelets 5 The Hartley Kernel On The Wavelet Analysis 

Based on the results obtained from Fourier-Like wavelets it is 
actually simple to define analytic wavelets. An analytic function 
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The Hartley transform kernel, or Hartley kernel, is defined by the 
“cosine and sine” function: cas(f) = cos(f) + sin(f). Recalling 



























































that Jf{sin }(/) = cos (t), the Hartley kernel can also be written 
as cas (t) = cos (t) —J4?{cos}(t) or cas (t) = sin(t) + J4?{sin}(t). 

This simple remark motivates the definition of Hartley-Like 
wavelets by taking the sum or the difference of a given real 
wavelet and its Hilbert transform. Let us define the Hartley ker¬ 
nel of a wavelet, or a Hartley-Like wavelet, < K{\i/(t)}, by: 

^{v(0} = -^-[v(0t^’{vH0]- ( 4 ) 

The following proposition proves that Tf{i f(t)} is also a 
wavelet and that the factor makes the Hartley kernel and its gen¬ 
erating wavelet to have same energy and admissibility coeffi¬ 
cient. Additionally, Proposition 6 shows that TC{ !//■(/)} has same 
null moments than their generating wavelet,! 

Proposition 5 Ify(t) is a real wavelet and its Hilbert 

transform, then < K{y/(t)} is also a real wavelet with same energy 
and admissibility coefficient of its generating wavelet, gff). 

Proof: If yf(t) is a real wavelet and f(t)} its Hilbert trans¬ 
form, then \jf(t) and J4?{\j/(t)} belong to L 2 (M) and hold the 
admissibility condition. The energy and admissibility coefficient 
of TC{!//■(/)} are given by: 

A simple manipulation yields to: 

1 r°° 

= ^Jjv(t)\ 2 +\^wm\ 2 dt, 

cmvKt)]= 

From Proposition[l] it follows that Ji{\j/(t)} also stays on L 2 (M) 
and it has same energy of \\f{t). Additionally, Tf{i f(t)} has also 
same admissibility coefficient than its generating wavelet, !//(/). 

□ 

Proposition 6 Let \\f{t) be a wavelet with N vanishing moments, 
then Tf{ has also N vanishing moments. 

Proof: From Proposition [2] it follows that the nth moment of 
< H{\f/(t)} is given by: 

L 1 n J (V2)"-(-2 -n-jY £o=q ’ 

which can also be written as 

M n [Ji{y}{t)\ = M n [yf(t)\ ±M n • 

Then, has also N null moments. 

In the frequency domain, Hartley-Like wavelets have the mag¬ 
nitude response of the generating wavelet multiplied by a scalar 
factor. Additionally, they impose a ±7r/4-shift on the phase re¬ 
sponse of the generating wavelet. 

Figure [2] shows some continuous-time wavelets and their cor¬ 
responding Hartley kernels using the addition operator in ([4]). 



Figure 3: A signal composed by two sinusoidal functions. Scale 
of the horizontal axis is 100 unit (time) per division, ranging 
from 0 to 1000. Scale of the vertical axis is 0.5 units (amplitude) 
per division, ranging from —2 to 2. 

6 Some Examples of Signal Analysis Using 
Fourier-Like and Hartley-Like Wavelets 

The wavelets proposed in this paper were simulated using the 
MATLAB Wavelet Toolbox © Standard sample signals were 
analyzed to illustrate the behavior of the proposed wavelets. 
Consider the wavelet transform coefficients C a ^ given by: 

1 r°° {t — b\ 

where a > 0 and b are real scale and translation scalars, respec¬ 
tively. 

6.1 Applying the Hartley-Like Wavelet Analysis 

Figure [3] shows a signal composed by two unitary-amplitude 
senoidal functions, 5 Hz and 9 Hz, and Figure [4] shows an 8- 
level wavelet analysis (wavelet transform) of that signal using 
the Morlet wavelet, its Hilbert transform, and its Hartley kernel. 

The scales where Hartley-Like Morlet wavelet may perform 
better signal analysis than the simple Morlet wavelet can be ob¬ 
served through the scales versus time charts. 

From Figure [4| significant differences on the wavelet trans¬ 
form at level 1 can be seen, when using an even, an odd or an 
asymmetrical wavelet. It follows that, despite the time shifting 
impelled by the wavelet transform, asymmetrical wavelets can 
retrieve more information from a signal regardless its symmetry. 
As the kind of symmetry of the analyzing signal is not a priori 
known, the use of Hartley-Like wavelets on the continuous-time 
wavelet transform may achieve better results. 

Figure [5] shows the signal derived from the lst-level wavelet 
transform of the later signal using this approach. Some improve¬ 
ments can be achieved using an odd or an asymmetrical wavelet 
to analyze an odd signal, than using an even wavelet. 

6.2 Applying the Analytic (or Fourier-Like) 
Wavelet Analysis 

Hilbert transform has already been utilized in wavelet-based sig¬ 
nal analysis {5}. In this section we carry a preliminary investiga- 
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(c) 




(d) 


(e) 


(f) 


Figure 2: Continuous-time wavelets and their Hartley kernels: (a) Morlet; (b) Meyer; (c) Mexican Hat; (d) Gaussian-1; (e) Gaussian- 
2; and (f) Gaussian-3. Scale of the horizontal axis is one unit (time) per division. Scale of the vertical axis is 0.2 units (amplitude) 
per division, ranging from — 1 to 1. 



(a) (b) (c) 


Figure 4: An 8-level wavelet analysis of later signal using: (a) the Morlet wavelet; (b) its Hilbert transform, and (c) its Hartley 
kernel. Scale of the horizontal axis is 100 unit (time or space b) per division, ranging from 0 to 1000. Scale of the vertical axis is 
one units (scale a) per division, ranging from 1 to 8. 
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Figure 5: The lst-level wavelet transform of a combined signal 
using the Morlet wavelet, its Hilbert transform and its Hartley 
kernel. Scale of the horizontal axis is 100 unit (time) per divi¬ 
sion, ranging from 0 to 1000. Scale of the vertical axis is 0.005 
units (amplitude) per division, ranging from —0.02 to 0.02. 



Figure 6: A standard frequency breakdown signal. Scale of the 
horizontal axis is 100 unit (time) per division, ranging from 0 
to 1000. Scale of the vertical axis is 0.2 units (amplitude) per 
division, ranging from — 1 to 1. 

tion about analytic wavelet analysis. Figure [6] shows a standard 
frequency breakdown signal and Figure [ 7 ] presents a 32-level 
wavelet analysis, assessed by the analytic Mexican Hat wavelet. 

From Figure [ 7 ] it can be observed that the Mexican Hat 
wavelet can identify the presence of both frequencies as well the 
time when occurs the frequency changing, which can be better 
accurate at lowest scale values. The modulus of that analytic 
wavelet analysis shows that the high frequency signal can be 
viewed at lowest scale values and the low frequency signal at 
higher scale values. 

As an example, Figures [8] and [9] show, respectively, the 5th 
and 25th-level wavelet transform, considering wavelets derived 
from the Mexican Hat wavelet. It can be seen the low and high 
frequency signals, when using real wavelets. 

Figure [T0| shows the normalized modulus of the 5th and 25th- 
level wavelet analysis, when using the analytic, or Fourier-Like, 
Mexican Hat wavelet. This feature shows that when analyzing 
individually scales, the time interval when occurs different fre¬ 
quency contents can be estimated. 



Figure 8: The continuous-time wavelet transform of a standard 
frequency breakdown signal using the Mexican Hat wavelet, its 
Hilbert transform, its Hartley kernel, and its Fourier kernel as 
mother wavelets and assuming 5 as the scale parameter. The 
modulus is used for the Fourier-Like wavelet analysis. Scale of 
the horizontal axis is 100 units (time) per division, ranging from 
0 to 1000. Scale of the vertical axis is one unit (amplitude) per 
division, ranging from —4 to 4. 



Figure 9: The continuous-time wavelet transform of a standard 
frequency breakdown signal using the Mexican Hat wavelet, its 
Hilbert transform, its Hartley kernel, and its Fourier kernel as 
mother wavelet and assuming 25 as the scale parameter. The 
modulus is used for the Fourier-Like wavelet analysis. Scale of 
the horizontal axis is 100 units (time) per division, ranging from 
0 to 1000. Scale of the vertical axis is one unit (amplitude) per 
division, ranging from —5 to 5. 
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(a) (b) (c) (d) 


Figure 7: A 32-level wavelet analysis of a standard frequency breakdown signal using the Analytic Mexican Hat wavelet: (a) real 
part of the wavelet coefficients; (b) imaginary part of the wavelet coefficients; (c) magnitude of the wavelets coefficients; and 
(d) phase of the wavelet coefficients. Scale of the horizontal axis is 100 units (time or scale b) per division, ranging from 0 to 100. 
Scale of the vertical axis is one unit (scale) per division, ranging from 1 to 31. 



Figure 10: Normalized modulus of the 5th and 25th-level 
wavelet coefficients of a standard frequency breakdown when us¬ 
ing the analytic, or Fourier-Like, Mexican Hat wavelet. Scale of 
the horizontal axis is 100 units (time) per division, ranging from 
0 to 1000. Scale of the vertical axis is 0.2 unit (amplitude) per 
division, ranging from 0 to 1. 

7 Conclusions 

New wavelet functions have been introduced that might improve 
the continuous-time signal analysis. Anti-symmetrical wavelets 
can be designed throughout the Hilbert transform of a sym¬ 
metrical wavelet, and vice-versa. Together, they are called a 


Hilbert transform pair of wavelets. Fourier-Like and Hartley- 
Like wavelets have been derived from the Fourier and Hart¬ 
ley transform kernels, which have been written, on the basis of 
Hilbert transform. 

The example cases illustrate that, despite the time shifting 
impelled by the wavelet transform, asymmetrical wavelets can 
grasp more features from a signal having no particular symme¬ 
try. Additionally, Fourier-Like and analytic wavelets have poten¬ 
tial applications for disturbances and frequency detection. 
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